Importance of the Electrostatic Correlations in Surface Tension of Hydrated Reline Deep Eutectic Solvent from Combined Experiments and Molecular Dynamics Simulations.

In this study, the surface tension and the structure of hydrated reline are investigated by using diverse methods. Initially, the surface tension displays a nonlinear pattern as water content increases, decreasing until reaching 45 wt %, then gradually matching that of pure water. This fluctuation is associated with strong electrostatic correlations present in pure reline, which decrease as more water is added. Changes in surface tension reflect a shift from charge layering in pure reline to an increased interfacial hydrogen bonding as the water content rises. This shift causes the segregation of urea molecules into the bulk phase and a gradual anchoring of water molecules to the air-reline interface. An interesting observation is the antisurfactant effect, where heightened interfacial anchoring results in an unexpected increase in real contribution of surface tension. This, along with weakened electrostatic correlations beyond 45 wt % due to reinforced interfacial hydrogen bonding, contributes to the complex behavior of surface tension observed in this study.

Microscopic insight into the shaping of MOFs and its impact on CO2 capture performance.

The traditional synthesis method produces microcrystalline powdered MOFs, which prevents direct implementation in real-world applications which demand strict control of shape, morphology and physical properties. Therefore, shaping of MOFs via the use of binders is of paramount interest for their practical use in gas adsorption/separation, catalysis, sensors, etc. However, so far, the binders have been mostly selected by trial-and-error without anticipating the adhesion between the MOF and binder components to ensure the processability of homogeneous and mechanically stable shaped MOFs and the impact of the shaping on the intrinsic properties of the MOFs has been overlooked. Herein, we deliver a first systematic multiscale computational exploration of MOF/binder composites by selecting CALF-20, a prototypical MOF for real application in the field of CO2 capture, and a series of binders that cover a rather broad spectrum of properties in terms of rigidity/flexibility, porosity, and chemical functionality. The adhesion between the two components and hence the effectiveness of the shaping as well as the impact of the overall porosity of the CALF-20/binder on the CO2/N2 selectivity, CO2 sorption capacity and kinetics was analyzed. Shaping of CALF-20 by carboxymethyl cellulose was predicted to enable a fair compromise between excellent adhesion between the two components, whilst maintaining high CO2/N2 selectivity, large CO2 uptake and CO2 transport as fast as in the CALF-20. This multiscale computational tool paves the way towards the selection of an appropriate binder to achieve an optimum shaping of a given MOF in terms of processability whilst maintaining its high level of performance.

Adsorption and dynamics of linear and mono-branched hexane isomers in MIL-140 metal-organic frameworks.

Recent breakthrough experiments revealed the iso-reticular Zr-MOFs, MIL-140B and MIL-140C, as promising sorbents for the separation of C6 isomers. Interestingly while the ultra-small pore MIL-140B exhibited hexane isomer sorption hierarchy according to the normal boiling point order (n-C6 > 3MP (3-methyl pentane)), an uncommon shift in the elution order was observed in the larger pore MIL-140C. It was only speculated that the flexibility of the MOFs might be the origin of this intriguing behavior. Herein, flexible force field hybrid osmotic Monte Carlo combined with molecular dynamics simulations were carried out to unravel the microscopic mechanism of the adsorption and dynamics of both C6 isomers in MIL140B and MIL140C. Thermodynamically preferred adsorption of n-C6 over 3MP was predicted for MIL-140B and to a slightly less extent for MIL-140C. Interestingly while the mobility of n-C6 was found to remain higher than that of 3MP in the whole range of loading for MIL-140B, 3MP becomes more mobile than n-C6 at saturation in MIL-140C. This suggests that this kinetics order is most probably the origin of the inversion of the elution order observed experimentally for MIL-140C. The translational and rotational dynamics of the two guests in MIL-140B and MIL-140C was further understood in-depth.

Molecular Dynamics Simulation Study of Organic Solvents Confined in PIM-1 and P84 Polyimide Membranes.

Organic solvent nanofiltration (OSN) has recently proved to be a promising separation process thanks to the development of membrane materials with suitable resistance toward organic solvents. Among those materials, P84 polyimide membranes are currently the most used in OSN while PIM-1 membranes have recently attracted attention due to their high permeance in apolar solvents and alcohols. Both P84 and PIM-1 membranes have nanosized free volumes, and their separation performance is finely connected to polymer/solvent interactions. Consequently, modeling OSN membranes at the molecular scale is highly desirable in order to rationalize experimental observations and gain a deeper insight into the molecular mechanisms ruling solvent and solute permeation. A prerequisite for understanding solvent transport through OSN membranes is therefore to characterize the membrane/solvent interactions at the molecular level. For that purpose, we carried out molecular simulations of three different solvents, acetone, methanol, and toluene in contact with P84 and PIM-1 membranes. The solvent uptake by both membranes was found to be correlated to the degree of confinement of the solvent, the polymer swelling ability and polymer/solvent interactions. The translational dynamics of the solvent molecules in the PIM-1 membrane was found to be correlated with the solvent viscosity due to the relatively large pores of this membrane. That was not the case with the P84 membrane, which has a much denser structure than the PIM-1 membrane and for which it was observed that the translational dynamics of the confined solvent molecules was directly correlated to the affinity between the P84 polymer and the solvent.

Cutting-edge molecular modelling to unveil new microscopic insights into the guest-controlled flexibility of metal-organic frameworks.

Metal-organic frameworks are a class of porous solids that exhibit intriguing flexibility under stimuli, leading often to reversible giant structural changes upon guest adsorption. DUT-49(Cu) and MIL-53(Cr) are fascinating flexible MOFs owing to their guest-induced breathing and negative gas adsorption behaviors respectively. Molecular simulation is one of the most relevant tools to examine these phenomena at the atomistic scale and gain a unique understanding of the physics behind them. Although molecular dynamics and Monte Carlo simulations are widely used in the field of porous materials, these methods hardly consider the structural deformation of a soft material upon guest adsorption. In this work, a cutting-edge osmotic molecular dynamics approach is developed to consider simultaneously the fluid adsorption process and material flexibility. We demonstrate that this newly developed computational strategy offers a unique opportunity to gain unprecedented molecular insights into the flexibility of this class of materials.

Molecular Insight into the Slow Dynamics of C4 Hydrocarbons in the Zeolitic-Imidazole Framework (ZIF-8).

The family of zeolitic-imidazole framework (ZIF) materials is currently considered for the challenging separation of C4 hydrocarbons. However, yet, the microscopic diffusion mechanism for these hydrocarbons in these narrow gate porous materials remains elusive by conventional methods due to its very slow nature. Experimental (solid-state 2H nuclear magnetic resonance-NMR) and computational (molecular dynamics-MD) approaches were applied together to derive slow diffusional dynamics of n-butane and 1-butene in ZIF-8. By means of the 2H NMR technique, we evidenced the presence of two adsorption sites for the guests localized inside the cages of ZIF-8 and in the vicinity of the gates. We characterized the molecular mobility at each site and revealed that the translational intercage diffusion is realized by a slow directional motion associated with the gate-crossing. MD simulations provide an in-depth analysis of the diffusion and fully support the proposed dynamics picture for both n-butane and 1-butene. These calculations enable the derivation of the diffusivity and barriers for the long-range diffusion of both hydrocarbons in ZIF-8 and unraveled the microscopic diffusion mechanism implying intracage and intercage motions. We show that this NMR approach combined with modeling is a valuable tool to probe the molecular mobility for slow diffusing species in ordered cagelike porous frameworks.

Dynamics of water confined in mesopores with variable surface interaction.

We have investigated the dynamics of liquid water confined in mesostructured porous silica (MCM-41) and periodic mesoporous organosilicas (PMOs) by incoherent quasielastic neutron scattering experiments. The effect of tuning the water/surface interaction from hydrophilic to more hydrophobic on the water mobility, while keeping the pore size in the range 3.5 nm-4.1 nm, was assessed from the comparative study of three PMOs comprising different organic bridging units and the purely siliceous MCM-41 case. An extended dynamical range was achieved by combining time-of-flight (IN5B) and backscattering (IN16B) quasielastic neutron spectrometers providing complementary energy resolutions. Liquid water was studied at regularly spaced temperatures ranging from 300 K to 243 K. In all systems, the molecular dynamics could be described consistently by the combination of two independent motions resulting from fast local motion around the average molecule position and the confined translational jump diffusion of its center of mass. All the molecules performed local relaxations, whereas the translational motion of a fraction of molecules was frozen on the experimental timescale. This study provides a comprehensive microscopic view on the dynamics of liquid water confined in mesopores, with distinct surface chemistries, in terms of non-mobile/mobile fraction, self-diffusion coefficient, residence time, confining radius, local relaxation time, and their temperature dependence. Importantly, it demonstrates that the strength of the water/surface interaction determines the long-time tail of the dynamics, which we attributed to the translational diffusion of interfacial molecules, while the water dynamics in the pore center is barely affected by the interface hydrophilicity.

Associated molecular liquids at the graphene monolayer interface.

We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.

Tuning the hexane isomer separation performances of Zeolitic Imidazole Framework-8 using mechanical pressure.

Hybrid osmotic Monte Carlo simulations were performed to anticipate the tunability of the separation performance of the flexible Zeolitic Imidazole Framework-8 (ZIF-8) via the application of an external mechanical pressure. This synergistic combination of mechanical control of the pore aperture/cage dimension and guest adsorption was applied to the challenging hexane isomers separation processes of vital importance in the field of petrochemical industry. The application of a mechanical pressure above 1 GPa was predicted to boost the linear hexane/2-methylpentane and 2-methylpentane/2,3-dimethylbutane selectivity by 40% and 17%, respectively, as compared to the pristine ZIF-8. We further unraveled the microscopic origin of this optimized performance with an in-depth analysis of the critical interplay between the structural changes of the ZIF-8 framework and the conformational rearrangements of C6 isomers under mechanical pressure.

Molecular Origin of the Prepeak in the Structure Factor of Alcohols.

Prepeak in the structure factor of alcohols is known for a half century and was attributed to one of two mechanisms (i) self-assembly in aggregates and (ii) existence of spatial heterogeneity. Although both explnations are often argued the molecular origin is yet unclear. In this work, molecular dynamics simulation of neat alcohols and their mixtures in the presence of an apolar liquid in bulk and in confined phases is performed to unveil and to clarify the origin of the prepeak at the molecular scale. Unambiguously, we show that the existence of the prepeak is the result of the self-assembly in clusters leading to long-range correlations rather than the spatial heterogeneity. We also establish that the confinement of neat liquids at the nanoscale does not erase the clustering and the prepeak but strongly reduce the spatial heterogeneity. Regarding the binary alcohol/toluene mixtures, we highlight the possibility to erase the clustering and the spatial heterogeneity from nanoconfinement inducing the formation of a core-shell structure. By tuning the interfacial chemistry and pore size, we shed light on the possibility to control the spatial heterogeneity, the self-assembly, and the microphase separation.

Force-Field Simulations of a Hydrated Lanreotide-Based Derivative: Hydration, Dynamics, and Numerical Evidence of Self-Assembly in Dimers.

Recently, self-organization of the cyclic octapeptide lanreotide and lanreotide-based derivatives in a nanotube to from a dimer structure has been experimentally evidenced. While the nature of the interactions between both monomers has been strongly investigated no molecular details of the hydration of the monomer and the formation of the dimer have been provided. Using molecular dynamics simulations, this work focuses on the structure, hydration, and dynamics of water and an analog of lanreotide. To do so, several models of monomers based on different schemes of partial charges and electrostatic interaction calculations are considered. By comparison with the experiments, we show that the model based on the combination of the AMBER force-field, CHELPG charge calculation, Ewald sum is the most relevant. Additionally, by mapping the interfacial hydration of the lanreotide monomer we evidence a heterogeneous surface in terms of hydrophilicity involving heterogeneous hydration. Furthermore, we show a slowdown in the translational dynamics of water molecules located close to the lanreotide surface. We also provide the molecular details of the self-assembly in the dimer in terms of structure, hydration, and energy.

Dynamic Heterogeneities in Liquid Mixtures Confined in Nanopores.

Binary liquid mixtures can exhibit nanosegregation, albeit being fully miscible and homogeneous at the macroscopic scale. This tendency can be amplified by geometrical nanoconfinement, leading to remarkable properties. This work investigates the molecular dynamics of tert-butanol (TBA)-toluene (TOL) mixtures confined in silica nanochannels by quasielastic neutron scattering and molecular dynamics simulation. It reveals a decoupling of the molecular motion of each constituent of the binary liquid, which can be followed independently by selective isotopic H/D labeling. We argue that this behavior is the signature of spatially segregated dynamic heterogeneities, which are due to the recently established core-shell nanophase separation induced by mesoporous confinement.

Tailoring the separation properties of flexible metal-organic frameworks using mechanical pressure.

Metal-organic frameworks are widely considered for the separation of chemical mixtures due to their adjustable physical and chemical properties. However, while much effort is currently devoted to developing new adsorbents for a given separation, an ideal scenario would involve a single adsorbent for multiple separations. Porous materials exhibiting framework flexibility offer unique opportunities to tune these properties since the pore size and shape can be controlled by the application of external stimuli. Here, we establish a proof-of-concept for the molecular sieving separation of species with similar sizes (CO2/N2 and CO2/CH4), via precise mechanical control of the pore size aperture in a flexible metal-organic framework. Besides its infinite selectivity for the considered gas mixtures, this material shows excellent regeneration capability when releasing the external mechanical constraint. This strategy, combining an external stimulus applied to a structurally compliant adsorbent, offers a promising avenue for addressing some of the most challenging gas separations.

Anomalous dynamics of water at the octopeptide lanreotide surface.

This work reports the study of water dynamics close to the cyclic octapeptide lanreotide from atomistic simulations of hydrated lanreotide, a cyclic octapeptide. Calculation of the hydrogen bonds between water molecules allows mapping of the hydrophilic regions of lanreotide. Whereas a super-diffusivity of the interfacial water molecules is established, a slowdown in rotational dynamics is observed, involving a decoupling between both processes. Acceleration in translation dynamics is connected to the hopping process between hydrophilic zones. Microscopically, this is correlated with the weakness of the interfacial hydrogen bonding network due to a hydrophobic interface at the origin of the interfacial sliding of water molecules. Heterogeneous rotational dynamics of water molecules close the lanreotide surface is evidenced and connected to heterogeneous hydration.

Anomalous Dynamics of a Nanoconfined Gas in a Soft Metal-Organics Framework.

Dynamics of confined molecules within porous materials is equally important as local structural order, and it is necessary to quantify it and to reveal the microscopic mechanisms ruling it for better control of adsorption applications. In this study, molecular dynamics simulations were carried out to investigate the translational and the rotational dynamics of methanol trapped into the flexible NH2-MIL-53(Al) metal-organics framework (MOF). Indeed, atomistic simulation is nowadays a relevant tool to explore matter at the nanoscale. Very recently it has been shown that the NH2-MIL-53(Al) MOF material was capable to undergo a reversible structural transition (breathing phenomenon) by combining adsorption and thermal stimuli. This flexibility can drastically affect the dynamics of confined molecules and therefore the successful conduct of adsorption applications such as gas storage and separation. Rotational and translational dynamics of confined methanol through nanoporous flexible NH2-MIL-53(Al) MOF were then deeply investigated by exploring a broad range of dynamical properties to extract the molecular mechanisms ruling them. This study allowed us to shed light on the interplay of dynamics of confined fluids and flexibility of porous material and to highlight the physical insights in diffusion mechanisms of confined molecules. Anomalous translational diffusion was evidenced due to a dynamical heterogeneity caused by a combination of a localized dynamics at the subnanometric scale and translational jumps between nanodomains in a zigzag scheme between the hydroxide group of the NH2-MIL-53(Al). Actually, the non-Fickian dynamics of methanol is the result of the specific host-guest interactions and the MOF flexibility involving the pore opening. Eventually, decoupling between both rotational and translational dynamics related to breaking in the Stokes-Einstein relation was highlighted.

Calculation of the interfacial tension of the graphene-water interaction by molecular simulations.

We report the calculation of the solid-liquid interface tension of the graphene-water interaction by using molecular simulations. Local profiles of the interfacial tension are given through the mechanical and thermodynamic definitions. The dependence of the interfacial tension on the graphene area is investigated by applying both reaction field and Ewald summation techniques. The structure of the interfacial region close to the graphene sheet is analyzed through the profiles of the density and hydrogen bond number and the orientation of the water molecules. We complete this study by plotting the profiles of the components of the pressure tensor calculated by the Ewald summation and reaction field methods. We also investigate the case of a reaction field version consisting in applying a damped shifted force in the case of the calculation of the pressure components.

Thermal and Guest-Assisted Structural Transition in the NH2-MIL-53(Al) Metal Organic Framework: A Molecular Dynamics Simulation Investigation.

Reversible structural transition between the Large (LP) and Narrow Pore (NP) forms (breathing phenomena) of the MIL-53(X, X = Al, Cr, Fe, Ga) Metal Organic Framework (MOF) is probably one of the most amazing physical properties of this class of soft-porous materials. Whereas great attention has been paid to the elucidation of the physical mechanism ruling this reversible transition, the effect of the functionalization on the flexibility has been less explored. Among functionalized MIL-53(Al) materials, the case of NH2-MIL-53(Al) is undoubtedly a very intriguing structural transition rarely observed, and the steadier phase corresponds to the narrow pore form. In this work, the flexibility of the NH2-MIL-53(Al) metal organic framework was investigated by means of molecular dynamics simulations. Guest (methanol) and thermal breathing of the NH2-MIL-53(Al) was thus explored. We show that it is possible to trigger a reversible transition between NP and LP forms upon adsorption, and we highlight the existence of stable intermediate forms and a very large pore phase. Furthermore, the NP form is found thermodynamically stable from 240 to 400 K, which is the result of strong intramolecular hydrogen bonds.

How Does the Surface Tension Depend on the Surface Area with Coarse-Grained Models?

We propose to investigate the size-effects on the surface tension calculated with coarse-grained (CG) models. We investigate different liquid-vapor (LV) and liquid-liquid (LL) interfaces with the MARTINI force field and original CG models designed for the dissipative particle dynamics (DPD) and multibody particle dynamics (MDPD) simulations. We also test a realistic CG potential developed for the DPD method to investigate the LV interface of n-pentane. Concerning the MARTINI force field, we observe a weak oscillatory effect of the interfacial tension with the surface area for the LV interfaces of n-octane and water. This weak dependence of the surface tension with the box dimension is also observed in the LL interface of n-octane-water (MARTINI, DPD) and in the LV interface of water with the MDPD model.

High Water Flux with Ions Sieving in a Desalination 2D Sub-Nanoporous Boron Nitride Material.

Over the past decades, desalination by reverse osmosis (RO) membranes has attracted increasing attention. Although RO has proven its efficiency, it remains, however, relatively costly because of the use of high-pressure pumps and the low water permeability of conventional cross-linked polymer membranes. One route to improve the desalination performance consists of using membranes made from sub-nanoporous boron nitride (sNBN) monolayers. Indeed, by using molecular dynamics simulations, we report here that the water permeability of such sNBN membranes far exceeds that of conventional RO polymer membranes and is even higher than that of nanoporous graphene while the ion rejection remains close to 100%. At the same time, the molecular mechanism of water and ion transport through sNBN has been elucidated, with special attention paid to the impact of ions on water permeability through sNBN membranes.